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Start-up  Response  of  Fluid  Film  Lubricated  Cryogenic 

Turbopumps  (Preprint) 

Luis  San  Andres1 

Mechanical  Engineering  Department, 

Texas  A&M  University,  College  Station,  TX  77843 

Reusable  primary  power  cryogenic  turbopumps  (TPs)  implement  externally  pressurized 
fluid  film  bearings  to  support  the  expected  large  thrust  and  lateral  radial  loads.  Compact  - 
low  count  part  TPs  operate  super  critically  at  exceedingly  high  shaft  speeds  (180  krpm)  and 
with  large  pressure  differentials.  Hybrid  journal  bearings  (HJBs)  enable  smaller  and  lighter 
turbopumps  through  no  bearing  DN  life  limitation.  The  growth  of  an  "all-fluid-film- 
bearing"  technology  for  reusable  and  less  costly  (per  launch)  TPs  demands  the  development 
of  analytical  models  and  design  tools  accompanied  by  the  testing  of  components.  The  paper 
presents  a  computational  model  for  the  prediction  of  the  start-up  performance  of  a  flexible 
rotor  supported  on  hydrostatic  radial  bearings.  The  transient  response  of  rotor-bearing 
systems  is  of  importance  to  determine  safe  operation  and  reliable  dynamic  performance 
under  extreme  loading  conditions.  In  the  start-up  operation  of  a  cryogenic  TP,  the  fluid 
supply  and  discharge  pressures,  as  well  as  the  radial  loads  acting  on  the  bearings,  depend  on 
pump  rotor  speed.  The  designed  aerodynamic  performance  of  the  whole  turbopump 
determines  the  schedule  of  rotor  speed  ramp-up.  The  start-up  event  is  quite  short  in  nature, 
lasting  a  few  seconds  at  most.  The  bearing  reaction  forces  are  calculated  from  the  numerical 
solution  of  unsteady  bulk-flow  equations  including  fluid  inertia,  flow  turbulence,  variable 
fluid  properties  and  thermal  energy  transport.  The  equations  of  motion  for  the  rotor  and 
the  fluid  film  bearings  are  solved  numerically  with  local  linearization  at  each  integration 
time  step.  Predictions  for  the  transient  start-up  response  of  a  test  rotor  supported  on  water 
hydrostatic  bearings  are  presented.  The  numerical  results  evidence  the  effect  of  rotor  mass 
on  the  rotordynamic  stability  of  the  rotor-bearing  system. 
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Nomenclature 


=  Cd7td02/4.  Effective  orifice  area  [m2] 
=  Recess  circumferential  length  [m] 

=  Radial  clearance  function  [m] 

=  Fluid  specific  heat  [J/kg  •  °K] 

=  Damping  force  coefficients  [Ns/m] 

=  Orifice  discharge  coefficient 
=  2  •  R.  Bearing  diameter  [m] 

=  Orifice  diameter  [m] 
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au  =  0.001375 
bM  =  500,000;cm  =  10.000 
eM  =  1/3.00 


Turbulent  flow  friction  factors  at  journal  and  bearing  surfaces 
=  Journal  center  eccentricity  components  [m] 

=  Bearing  fluid  film  forces  along  {X,  Y}  axes  [N] 

=  c  +  X(t)  cos(9)  +  Y(t)  sin(O).  Film  thickness  [m] 

=  Recess  depth  [m] 

=  Heat  convection  coefficients  [J/kg  •  °K] 

=  Bearing  force  stiffness  coefficients  [N/m] 

=  Bearing  axial  length,  recess  axial  length  [m] 
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=  Bearing  inertia  force  coefficients  [kg] 

=  Fluid  pressure,  recess  pressure,  supply  pressure  [N/m2] 

=  Heat  flow  into  bearing  and  journal  (shaft)  surface  [J/kg] 

=  ( pQCRJ  p )  .  Nominal  circumferential  flow  Reynolds  number 

=  ( p/p)h-\[ (vx  -QR)2  +  vl] ;  (p/p)h^[V2x  +  V2J 

Flow  Reynolds  numbers  relative  to  journal  and  bearing  surfaces 
=  Surface  roughness  depths  at  journal  and  bearing  [m] 

=  Time  [s] 

=  Temperature,  supply  temperature  [°K] 

=  Bulk-flow  velocities  in  circumferential  (x)  and  axial  (z)  directions  [m/s] 

=  Recess  volume  including  supply  line  [m3] 

=  External  loads  applied  on  journal  [N] 

=  Fluid  swirl  ratio  at  recess  edges 

=  (l/p)(dp/dP).  Liquid  compressibility  coefficient  [m2/N] 

=  -(l/p)(dp/dT).  Liquid  volumetric  expansion  coefficient  [1/°K] 

=  x/R.  Circumferential  or  angular  coordinate. 

=  Zi  (kj  +  kb).  Turbulence  shear  factors  along  axial  and  circumferential  flow  directions 
=fj  ■ Rej,fB  ■  ReB  .  Turbulent  shear  parameters  at  journal  and  bearing  surfaces 
=  Fluid  density  [kg/m3],  viscosity  [Ns/m2] 

=  Empirical  recess-edge  entrance  loss  coefficients  in  circumferential  (upstream,  downstream) 
direction 

=  Empirical  recess-edge  entrance  loss  coefficients  in  axial  direction 
=  Rotational  speed  of  journal,  excitation  or  whirl  frequency  [1/s] 

=  In  direction  of  local  circumferential  and  axial  coordinates  in  plane  of  bearing 
=  Orifice 

=  Bearing  recesses  and  edges  (entrance) 

=  Upstream  and  downstream  of  recess 
=  Refer  to  bearing  and  journal  surfaces 


I.  Introduction 

Hydrostatic  bearings  (HBs)  derive  their  load  capacity  not  from  shear  flow  driven  effects  (hydrodynamic  wedge 
and  surface  sliding)  but  rather  from  the  combination  of  pressure  versus  flow  resistance  effects  through  a  feed 
restrictor  and  within  the  bearing  film  lands.  Hydrostatic  bearings  can  support  large  loads  without  journal  rotation 
and  provide  large  (accurate  and  controllable)  direct  stiffness  as  well  as  damping  coefficients.  The  hydrostatic 
stiffness  is  of  unique  importance  for  the  centering  of  high-precision  milling  machines,  gyroscopes,  large  arena 
movable  seating  areas,  telescope  bearings,  and  even  cryogenic  fluid  turbo  pumps  for  rocket  engines. 

The  importance  of  hybrid  (combination  hydrostatic  and  hydrodynamic)  journal  and  thrust  bearings  and  damping 
seal  bearings  as  radial  support  elements  in  cryogenic  turbomachinery  has  steadily  grown  over  the  past  few  years1'3. 
Compact  -  low  count  part  turbo  pumps  (TPs)  operate  sub  critically  at  exceedingly  high  shaft  speeds  (180  krpm)  with 
pressure  differentials  as  large  as  550  bar.  Advanced  primary  power  require  of  externally  pressurized  fluid  film 
bearings  to  support  the  expected  large  thrust  and  lateral  radial  loads.  Hybrid  journal  bearings  (HJB)s  enable  smaller 
and  lighter  turbopumps  through  no  bearing  DN  life  limitation  and  sub  critical  rotor  operation,  i.e.  at  speeds  below 
the  first  elastic  mode  of  the  rotor-bearing  system.  HJBs  offer  durability,  low  friction  and  wear,  accuracy  of 
positioning,  and  large  direct  stiffness  and  damping  force  coefficients.  These  features  enable  the  design  (and 
operation)  of  un-shrouded  impellers  with  a  significant  increase  in  the  turbopump  mechanical  efficiency.  The  growth 
of  an  "all-fluid-film-  bearing"  technology  for  advanced  and  less  costly  (per  launch  cost)  turbopumps  demands  the 
development  of  analytical  models  and  design  tools,  the  testing  of  components,  and  the  implementation  of  the 
technology. 

Hydrostatic  bearings  require  an  external  pressurized  supply  system  and  some  type  of  flow  restrictor.  Also, 
under  dynamic  motions,  hydrostatic  bearings  may  display  a  pneumatic  hammer  effect  due  to  fluid  compressibility. 
However,  and  most  importantly,  the  load  and  static  stiffness  of  a  hydrostatic  bearing  are  independent  of  fluid 
viscosity;  thus  making  this  bearing  type  very  attractive  for  cryogenic  liquid  turbopumps  or  low  viscosity  process 
fluid  pump  applications. 
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Note  that  for  the  cryogenic  fluid  application  as  well  as  others  handling  low  viscosity  liquids,  the  large  surface 
speeds  and  the  large  pressure  differential  determine  flow  conditions  with  high  levels  of  flow  turbulence  and  fluid 
inertia  effects.  Flow  turbulence  increases  the  lubricant  “effective”  viscosity,  thus  enhancing  the  load  capacity  due  to 
hydrodynamic  effects  and  increasing  the  bearing  energy  dissipation  characteristics,  i.e.  more  damping12. 
Computational  programs  based  on  the  Reynolds  equation  of  classical  lubrication,  i.e.  no  fluid  inertia,  are  ill-prepared 
to  render  adequate  predictions  of  hybrid  bearing  performance,  static  and  dynamic  force  coefficients. 

The  author  has  developed  the  most  comprehensive  computational  analyses  for  prediction  of  process  fluid  hybrid 
bearings,  radial  and  thrust.  The  analyses  address  to  the  most  important  theoretical  and  practical  issues  related  to  the 
operation  and  dynamic  performance  of  cryogenic  fluid  film  bearings,  i.e.  geometric  configuration,  operating 
conditions,  flow  turbulence,  fluid  inertia,  realistic  fluid  properties,  thermal  effects,  and  two-phase  flow  phenomena. 
References  1,  2  and  4-12  detail  the  computational  analyses  along  with  measurements  to  validate  and  calibrate  the 
predictive  codes.  San  Andres13  presents  a  comprehensive  review  of  the  static  and  dynamic  performance 
characteristics  of  annular  pressure  seals  and  hydrostatic  bearings  for  process  fluid  high  performance  pumps.  The 
major  advantages  and  disadvantages  of  hydrostatic  bearings  are  thoroughly  discussed  with  emphasis  on  remedies  or 
fixes  to  reduce  or  even  eliminate  the  potential  for  hydrodynamic  instability  and  pneumatic  hammer  instability  in 
actual  TP  applications. 

The  objective  of  the  present  analysis  is  to  advance  a  computational  model  for  the  prediction  of  the  start-up  (or 
shut  down)  performance  of  a  rotor  supported  on  hydrostatic  bearings.  In  the  start-up  operation  of  a  cryogenic 
turbopump,  the  fluid  supply  and  discharge  pressures,  as  well  as  the  radial  loads  acting  on  the  hydrostatic  bearings, 
depend  on  the  pump  rotor  speed.  The  schedule  of  rotor  speed  ramp-up  is  defined  by  the  designed  aerodynamic 
performance  of  the  whole  turbopump.  The  start-up  event  is  quite  short  in  nature,  lasting  a  few  seconds  at  most. 

The  transient  response  of  rotor-bearing  systems  is  of  importance  to  determine  safe  operation  and  dynamic 
performance  under  extreme  loading  conditions.  San  Andres14  already  advanced  a  simpler  model  for  the  transient 
response  of  a  point-mass  rotor  supported  in  turbulent  flow,  externally  pressurized  fluid  film  bearings.  The  equations 
of  motion  are  solved  numerically  with  local  linearization  at  each  integration  time  step.  The  bearing  reaction  forces 
are  calculated  from  the  numerical  solution  of  unsteady  bulk-flow  equations  including  fluid  inertia,  turbulence, 
variable  fluid  properties  and  thermal  energy  transport.  Examples  follow  for  the  transient  response  of  damper  seals 
and  hydrostatic  bearings  under  a  variety  of  external  loads.  The  nonlinear  model  needs  small  time  steps  with  large 
execution  times.  Note  that  the  analysis  in  Ref.  14  is  limited  to  predictions  of  transient  response  at  a  fixed  rotor 
speed.  Presently,  this  severe  shortcoming  is  removed. 

II.  Unsteady  bulk  flow  analysis  of  turbulent  flow  hydrostatic  bearings 

Figure  1  shows  the  geometry  of  a 
hybrid  (combination  hydrostatic 
/hydrodynamic)  journal  bearing.  A 
liquid  at  high  pressure  ( Ps )  is  supplied 
through  orifice  restrictors  and  impinges 
into  the  bearing  recesses  with  a  mean 
pressure  (PR).  The  fluid  injection  is 
typically  radial;  though  in  some  instances 
it  could  be  at  an  angle  opposing  shaft 
rotation2.  The  pressure  field  within  the 
recesses  is  determined  from  flow 
continuity  with  the  film  lands, 
momentum  exchange  at  the  orifice  plane 
and  a  viscous  rise  due  to  journal  rotation. 
At  the  recess  edges,  an  inertial  pressure 
drop  also  occurs  due  to  the  sudden 
transition  from  the  recess  of  depth  (hR) 
into  the  film  lands  of  thickness  ( h ).  Past 
the  recesses,  the  liquid  then  flows 
through  the  film  lands  and  the  pressure 


Figure  1.  Schematic  view  of  a  radial  hydrostatic/hydrodynamic 
journal  bearing 


2  Angled  injection  aids  to  reduce  the  development  of  circumferential  flow  speed  and  reduce,  even  eliminate,  the 
magnitude  of  cross-coupled  stiffness  coefficients11. 
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drops  to  the  discharge  value  (Pa). 


The  computational  model  considers  the  fully  developed  turbulent  bulk-flow  of  a  fluid  whose  material  properties 
depend  on  its  local  thermo  physical  state  of  pressure  and  temperature.  That  is,  the  fluid  density  (p ),  viscosity  (ju )  and 
heat  capacity  (CP)  are  functions  of  the  local  pressure  (P)  and  temperature  (7).  The  fluid  flows  within  the  film  lands 
of  film  thickness  (/?)  with  buk-flow  (film  averaged)  circumferential  and  axial  flow  velocities  denoted  by  Vx  and  V:, 
respectively.  The  general  bulk-flow  transport  equations  are1: 


d{phy/)  |  d{phVxy/)  |  d{phV:y/)_^ 


dt 

dx 

dz 

where 

Variable 
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Above,  (kx  k-J  denote  wall  shear  stress  turbulent  flow  coefficients.  The  wall  shear  stress  parameters  /<; = kx  =  Zi ( Kt  +  kb ) 
with  Kj=fj  Rej,  KB=fB  ReB,  and  the  Moody’s  friction  factors  (Jj.n)  depend  on  the  bearing  and  journal  surface 
conditions  and  the  flow  Reynolds  numbers,  Rej  and  ReB,  relative  to  the  rotating  and  stationary  surfaces, 
respectively13.  The  bearing  film  thickness  relates  the  rotor  motion  at  the  journal  location  to  the  bearing  fluid  flow, 

h=c(x)  +  Xftj  cos(<9)+  y(l]  sin(6»)  (3) 


and  its  time  derivative,  dh/dt.  Above,  c  is  the  bearing  clearance  at  the  journal  centered  position,  and  ( X ]  Y)  are  the 
instantaneous  journal  center  displacements.  The  journal  eccentricity  follows  from  e=(X2+Y2)112.  In  the  transport  of 
thermal  energy,  Eq.  (2d),  QBS  =QB  +  Qs  =  hB  (T-TB)+hs  (T-Ts)  is  the  heat  flow  through  the  bearing  and  shaft 
surfaces  at  temperatures  TB  and  Ts,  respectively  .  fir  is  the  coefficient  of  thermal  expansion,  (hB  ,  hs)  are  heat 
convection  coefficients,  and  /lrz>,  and  Azvf/  are  wall  shear  stress  differences  along  the  axial  and  circumferential 
directions,  respectively. 

In  most  time  transient  response  analyses,  the  time  scales  for  thermal  energy  transport  and  diffusion  are  much 
slower  than  those  for  mechanical  energy  (due  to  vibrations,  for  example).  Hence  a  quasi-state  energy  transport  is 


assumed,  with  the  temperature  remaining  constant  during  a  relatively  fast  transient  process,  i.e.  ST  /  dt  ~  0 

Figure  2  depicts  a  bearing  recess 

with  axial  length  (/)  and 

circumferential  extent  ( b ).  The  recess 

area  AR=l  x  b,  and  the  feed  orifice  has 

diameter  da  with  feed  volume  Vsuppir 

The  simplified  analysis  of  hydrostatic 

bearings  does  not  model  the  flow  field 

within  the  recess  since  these  are 

(typically)  deep  and  enclose  a  nearly 

stagnant  fluid  volume.  The  bulk-flow 

model  accounts  for  mass  flow 

„  ,  _  .  ..  ,  .  ...  .  .  .  .  .  continuity  with  the  film  lands  and 

Figure  2.  Schematic  view  of  approximate  pressure  field  in  a  pocket  ,  n  ^  c 

r  .  .  ,  , .  .  obtains  the  recess  pressures  ( PR )  from 

of  a  hydrostatic  bearing 
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an  orifice  flow  equation  whose  accuracy  depends  on  an  empirical  discharge  coefficient  (C,/).  CFD  results  and 
measurements  show  the  generation  of  hydrodynamic  pressures  within  the  pocket,  followed  by  sharp  inertial  pressure 
drops  at  the  recess  edges. 

The  continuity  equation  at  a  hydrostatic  recess  establishes  a  balance  among  the  mass  flow  through  the  feed 
orifice  (Mr),  the  flow  through  the  boundaries  of  the  recess  into  the  film  lands  (Mr),  and  the  accumulation  of  fluid 
mass  within  the  recess  volume,  VR=[AR  (h+hg)+V suppiy] .  That  is, 

(  Y/2  d 

MR=CdA0\£\P,-PR\  J  =MY+—{pVR) 

where  A0  =  Cj  ( jtc102/4)  is  the  effective  orifice  area,  and  ^  =  ^ph  V  -rj)dr  *s  outflow  from  the  pocket  into  the 

r 

bearing  film  lands.  The  circumferential  pressure  downstream  of  the  feed  orifice,  PR  ,  increases  as  in  a  Rayleigh  step 
bearing,  i.e. 


Pr  =  Pr+Pkx 


2  ( h  +  hR  )2 


Q  R 


-V, 


(5) 


Fluid  inertia  causes  a  sudden  pressure  drop  at  the  interface  between  a  recess  and  the  film  lands.  The  fluid  pressures, 
PR  ,  entering  into  the  film  lands  bounding  a  recess  are 


p«=K^P 


1- 


V;v 


Pe  ) 


\'l  +  llR  j 


v: 


(6) 


where  ( c)  represents  empirical  entrance  loss  coefficients  at  the  recess  edges,  axial  and  circumferential.  The  sudden 
pressure  drop  is  accounted  for  only  if  the  fluid  flow  effectively  enters  into  the  thin  film  lands. 

Recall  that  severe  subsynchronous  vibrations  at  rotational  speeds  above  a  certain  threshold  denote  a 
hydrodynamic  instability  on  rotor- fluid  film  bearing  systems  and  due  to  the  effect  of  journal  rotational  speed  on  the 
shear  flow  field.  This  condition  is  typical  of  fixed  geometry  bearings.  The  threshold  speed  corresponds  to  the  rotor 
speed  at  which  a  bearing  is  deprived  from  its  effective  damping  and  any  small  perturbation  from  an  equilibrium 
position  will  determine  unbounded  rotor  motions.  The  whirl  frequency  ratio  ( WFR )  denotes  the  ratio  between  the 
onset  whirl  frequency  (typically  the  system  first  critical  speed)  and  the  threshold  speed  of  instability.  Plain  journal 
bearings  show  a  WFR  equal  to  0.50  for  small  to  moderate  operating  eccentricities  (light  loads),  and  thus  instability 
onsets  at  rotational  speeds  equal  to  twice  the  system  first  critical  speed.  Measurements  in  hybrid  bearings  verify 
closely  the  prediction  of  WFR  =0.50,  see  Ref.  10.  In  some  circumstances  the  WFR  even  increases  above  0.50,  in 
particular  for  low  rotational  speeds  and  large  supply  pressures. 

San  Andres  and  Childs"  extend  the  bulk-flow  model  to  account  for  fluid  injection  at  an  angle  and  opposing 
shaft  rotation.  This  design  feature  retards  the  frill  development  of  the  circumferential  flow  velocity,  thus  reducing 
the  cross-coupled  stiffness  coefficients  which  prevent  the  operation  of  hybrid  bearings  at  large  rotational  speeds. 

ITT.  Numerical  solution  of  unsteady  bulk-flow  equations 

In  the  numerical  analysis,  an  implicit  scheme  is  implemented  on  finite  size  control  surfaces16  with  the  SIMPLEc 
algorithm17.  In  general,  the  discrete  algebraic  equation  for  each  variable  <p=  {V„  V:,  P,  T},  at  the  current  time  (/),  is 
of  the  form 

A*4>p=A*4>E+A*r4>w+A*4>s+AtiN+S*  +B*^St  (7) 

where  ’’  is  the  value  of  the  variable  one  time  step  before.  That  is,  at  each  time  step,  the  previous  flow  field  must 
be  known  fully;  in  particular  the  one  at  the  initial  time  when  the  solution  procedure  starts.  San  Andres18  details  the 
development  of  the  algebraic  equations  for  momentum,  mass  and  energy  transport  in  thin  film  bulk-flows.  Once  the 
solution  of  the  set  of  Eq.  (7)  is  obtained  at  time  t ,  fluid  film  bearing  reaction  forces  (Fx,  Fr)  are  calculated  from 
integration  of  the  pressure  field,  i.e. 
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RdOdz 


(8) 


-}  J  (p-p.) 


cos  0 
sin0 


The  bearing  reaction  forces  are  also  known  as  (nonlinear)  impedances  since  they  are  a  function  of  the 
instantaneous  journal  center  position  and  its  velocity  components,  i.e. 

Fx=fx(x,Y,X,  y)  Fy  =  fY  (X,  Y, X,  V)  (9) 

A  rotordynamics  model  predicting  the  transient  response  of  a  rotor  supported  on  bearings  needs  to  integrate  the 
(nonlinear)  bearing  impedances  at  each  time  step.  The  typical  equations  of  motion  are  of  the  form 

[M]  {ii}-[[C  ]  -  Q  [C]R  ]  \Ct)+  [k]{u}=  {Fext  (u,ii,t)}+  {Fb  (i u,u )}  (,()) 

where  O  is  the  rotor  speed,  and  [M\,  [  K\.  [C]  and  [G]  denote  the  system  mass,  stiffness,  damping  and  gyroscopic 
matrices.  The  vector  {u}  represents  rotor  displacements  (translations  and  rotations),  {F  t  (it ,  u,t )}  denotes  external 

forces  such  as  static  (gravity)  loads  and  those  due  to  mass  imbalance,  and  {Fb  ( u ,  it ) }  corresponds  to  the  nonlinear 
bearing  impedances,  for  example,  i.e.  Eq.  (8).  Note  that  the  solution  of  the  rotor-bearing  system  equations  of  motion, 
Eq.  (10),  is  linked  to  the  solution  of  the  bulk-flow  equations  for  each  bearing  support,  Eq.  (7). 

The  interface  of  the  current  predictive  hydrostatic  bearing  bulk-flow  model  to  a  comprehensive  (realistic) 
rotordynamics  model  is  being  completed.  At  this  time,  for  testing  of  the  code  features,  a  simple  point  mass  rotor¬ 
bearing  model  suffices.  The  equations  of  motion  for  the  simple  system  are 

MX  =  WX  +FX,  MY  =  Wy+Fy  (11) 

where  M  denotes  the  rotor  mass,  ( Wx,  Wy)  and  (FX,  Fy)  are  the  components  of  the  externally  applied  force  and  the 
bearing  reaction  forces,  respectively. 

The  numerical  solution  of  Eq.  (11)  implements  the  implicit  Wilson-0  method  as  given  by  Bathe19  and 
reproduced  in  Ref.  14  with  0=1.42.  At  each  time  step,  Eq.  (11)  is  an  algebraic  nonlinear  equation  whose  unique 
solution  is  determined  with  the  aid  of  linearized  stiffness  and  damping  coefficients. 

IV.  Predictions  of  transient  speed  response  of  a  simple  rotor  supported  on  water  lubricated 

hydrostatic  bearings 

A  unique  test  rig  has  been  constructed  at  Texas  A&M  University  for  purposes  of  recording  the  start  up  response 
of  a  rotor  supported  on  a  water  lubricated  hydrostatic  bearing.  The  test  apparatus  will  replicate  a  typical  rotor  speed 
start-up  in  a  cryogenic  TP  and  measure  the  operating  conditions  best  suited  for  early  rotor  lift  off.  Shallow  pockets 
of  small  area  are  preferred  in  high  performance  TPs  to  avoid  pneumatic  hammer  effects  when  handling  compressible 
liquids.  Details  of  the  specific  bearing  configuration  are  omitted  for  brevity.  In  the  experiments,  warm  water  will  be 
supplied  into  the  bearings  with  increasing  supply  pressures  tied  to  the  rotor  speed  ramp-up.  For  the  current 
predictions,  isothermal  flow  conditions  prevail  in  the  bearings  since  the  flow  rates  are  large  and  the  shaft  speed-up 
ramp  is  very  fast.  That  is,  in  the  analysis  the  fluid  temperature  remains  at  its  inlet  condition. 

Figure  3  shows  the  schedule  of  a  fast  transient  for  the  envisioned  test  rig.  The  figure  displays  the  rotor  speed, 
water  supply  pressure  and  external  radial  load  acting  on  the  bearing  versus  time.  For  times  less  than  0.015  sec,  the 
rotor  speed  increases  linearly  with  time  from  5  kipm  to  20  kipm.  Over  the  same  time  interval,  the  radial  load  also 
increases  linearly  from  100  N  to  300  N,  while  the  supply  pressure  has  a  quadratic  increase,  i.e.  proportional  to  rotor 
speed2.  For  times  larger  than  0.015  sec,  the  speed,  feed  pressure  and  load  remain  invariant.  Figure  4  depicts  the  feed 
pressure  and  radial  load  (Wx,  fFy=0)  versus  rotor  speed.  Over  the  time  span  from  0  s  to  0.027  s,  the  total  number  of 
shaft  revolutions  is  just  seven. 

In  the  transient  rotor  response  (nonlinear)  analysis,  a  time  step  At=33.3  micro-s  is  used,  corresponding  to  a 
sampling  rate  of  30  kHz.  For  the  lowest  and  top  shaft  speeds  of  5  krpm  and  20  krpm  (83  Hz  &  333  Hz),  the  number 
of  steps  per  shaft  revolution  equals  to  360  and  90,  respectively.  The  fast  speed  ramp  rate,  15  krpm/0.015  sec  (16.6 
kHz/sec)  demands  the  very  small  time  step  for  integration  of  Eq.  (11). 
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Figure  4:  Example  transient  supply  pressure  and 
static  load  vs  rotor  speed  for  water  HB  test  rig 


Figure  3.  Example  transient  rotor  speed, 
supply  pressure  and  static  load  vs.  time  for 
water  HB  test  rig 

Numerical  predictions  of  hydrostatic  bearing  static  performance  were  obtained  prior  to  conducting  the 
integration  of  the  rotor-bearing  system  transient  response  during  the  ramp  up  in  shaft  speed.  Table  1  depicts  the 
values  of  speed,  supply  pressure  and  static  load  acting  on  the  bearing,  and  the  predicted  rotor  static  position  (Xs,  Ys) 
mass  flow  and  power  loss,  equivalent  stiffness,  whirl  frequency  ratio  and  critical  mass  at  the  noted  shaft  speeds. 
Note  that  the  whirl  frequency  ratio  ( WFR )  is  >~  0.50  for  all  speeds,  denoting  the  potential  of  the  bearings  to  induce 
a  hydrodynamic  instability.  The  critical  mass  parameter  decreases  rapidly  as  the  rotor  speed  increases.  The  threshold 
speeds  of  instability  are  also  noted  for  two  rotor  masses  ( M)  equal  to  1  and  10  kg.  For  the  largest  rotor  mass,  a 
threshold  speed  of  20.36  krpm  is  very  close  to  the  top  shaft  speed  (20  krpm)3. 

As  an  initial  condition  for  the  numerical  integration  of  the  equations  of  motion,  the  flow  field  corresponding  to 
the  journal  static  position  at  the  lowest  shaft  speed  (5  krpm)  was  used,  i.e.  A7c=0.38,  Y/c= 0.15.  Figures  5and  6 
depict  the  predicted  dimensionless  (with  respect  to  bearing  clearance)  journal  eccentricity  and  its  components  (X,Y) 
versus  shaft  speed  and  versus  time,  respectively.  The  graphs  include  the  results  for  rotor  mass  equal  to  1  kg  and  10 
kg.  For  the  smallest  rotor  mass,  M=  1  kg,  the  dynamic  response  predictions  agree  well  with  the  static  load 
calculations;  and  most  importantly  show  a  stable  rotor-bearing  system.  On  the  other  hand,  for  the  largest  rotor  mass, 
M=10  kg,  the  transient  response  quickly  evolves  into  an  oscillatory  pattern  with  increasing  amplitudes  of  motion, 
typical  of  the  onset  of  a  hydrodynamic  instability.  The  frequency  of  whirl  is  subsynchronous,  tracking  the  rotor 
speed  at  roughly  55%  WFR. 


The  predictions  for  critical  mass  and  threshold  speed  of  instability  do  account  for  fluid  inertia  effects  in  the 
bearings.  This  effect  amounts  to  27%  of  the  static  direct  stiffness  at  the  top  shaft  speed. 
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(a)  rotor  mass,  M=  1  kg 


-  Transient 

□  □  Static 


(b)  rotor  mass,  M=  10  kg 

Static  &  Dynamic  Response 


shaft  speed  (RPM) 


Locus  of  journal  center 


nnD  Static 


Figure  5.  Predicted  dimensionless  journal  (rotor)  eccentricity  versus  shaft  speed  (a)  M=1  kg,  (b)  M=10  kg. 
Static  load  journal  eccentricity  shown  as  square  symbol. 


....  y/c 
-  e/c 

Figure  6a:  Predicted  dimensionless  journal  (rotor)  eccentricity  components  versus  time.  Rotor  mass  M=  1  kg 
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---•  Y/c 
-  e/c 

Figure  6b.  Predicted  dimensionless  journal  (rotor)  eccentricity  components  versus  time.  Rotor  mass  M=  10  kg 


Table  1:  Static  load  and  feed  pressure  and  predicted  rotor  (journal)  static  position,  flow  parameters,  critical 

mass  and  threshold  speed  of  instability 


X/c 

Y/c 

Mass  Flow 

Speed 

P  Supply 

P  Exit 

Load-X 

Load-Y 

rpm 

bars 

bars 

N 

N 

H 

[-1 

kg/s 

5000 

5.0 

1.0 

100 

0 

0.38 

0.15 

0.16 

10000 

6.7 

1.0 

167 

0 

0.37 

0.23 

0.18 

15000 

11.7 

1.0 

233 

0 

0.27 

0.18 

0.26 

20000 

20.0 

1.0 

300 

0 

0.19 

0.14 

0.36 

10  kg 

1  kg 

Speed 

Power  Loss 

Keq 

WFR 

Critical  mass 

Threshold 

speed 

Threshold 

speed 

rpm 

kW 

N/m 

- 

Kg 

rpm 

rpm 

5000 

0.01 

3.53E+06 

0.50 

48.8 

11041 

34915 

10000 

0.06 

4.80E+06 

0.54 

12.4 

11129 

35193 

15000 

0.17 

9.44E+06 

0.53 

10.8 

15576 

49255 

20000 

0.37 

1.66E+07 

0.53 

10.4 

20356 

64371 

Speed 

Kxx 

Kxy 

Kyx 

Kyy 

Cxx 

Cxy 

Cyx 

cyy 

Mxx 

Mxy 

Myx 

Myy 

rpm 

N/m 

N/m 

N/m 

N/m 

N-s/m 

N-s/m 

N-s/m 

N-s/m 

kg 

kg 

kg 

kg 

5000 

3.59E+06 

1.32E+06 

-7.54E+05 

3.16E+06 

3747 

484 

-621 

3783 

0.9 

-0.2 

0.0 

1.0 

10000 

4.13E+06 

2.70E+06 

-2.67E+06 

4.35E+06 

4830 

1101 

-893 

4723 

0.8 

-0.1 

0.0 

0.8 

15000 

8.12E+06 

5.42E+06 

-5.38E+06 

8.27E+06 

6577 

1557 

-1448 

6494 

0.8 

-0.1 

0.1 

0.9 

20000 

1 .43E+07 

9.65E+06 

-9.53E+06 

1 .43E+07 

8681 

2105 

-2037 

8613 

0.9 

-0.1 

0.1 

0.9 

Figures  7  and  8  present  the  applied  external  load  (Wx)  and  predicted  FIB  reaction  force,  IT7!,  versus  shaft  speed 
and  versus  time,  respectively.  For  the  lowest  rotor  mass,  M=\  kg,  inertial  effects  are  negligible;  and  hence  the 
motion  is  quasi-static,  i.e.  for  all  times,  WX+FA= 0  and  Fy~0.  On  the  other  hand,  for  the  large  rotor  mass,  M=10  kg, 
the  bearing  reaction  force  shows  a  superposition  of  the  static  force  (Wx)  and  a  periodic  force  whose  magnitude 
increases  steadily  as  the  shaft  speed  increases.  The  frequency  of  the  dynamic  force  component  is  subsynchronous. 

In  summary,  the  numerical  predictions  obtained  for  the  ramp  in  shaft  speed  with  varying  feed  pressure  into  the 
hydrostatic  bearing  and  increasing  radial  loads  acting  on  the  rotor  are  in  agreement  with  prior  knowledge.  Most 
importantly,  the  predictions  evidence  the  importance  of  rotor  mass  on  the  response  of  the  system.  That  is,  for  a  small 
rotor  mass,  “rotor  inertia  forces”  are  negligible,  and  the  hydrostatic  bearing  reacts  quasi-statically  to  the  applied 
external  force  on  the  rotor.  On  the  other  hand,  for  a  larger  rotor  mass,  the  system  becomes  unstable  with  a  typical 
-50%  whirl  frequency.  This  phenomenon  is  predicted  accurately  from  a  linear  stability  analysis  of  the  rotor-bearing 
system  using  the  bearing  force  coefficients  (stiffness,  damping  and  inertia);  albeit  the  amplitude  of  motion  can  only 
be  determined  from  the  solution  of  the  nonlinear  equations  of  motion. 
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(a)  rotor  mass,  M=  1  kg 


Static  &  Dynamic  Response 


□  □  Static 


(b)  rotor  mass,  M=  10  kg 


shaft  speed  (RPM) 

-  Transient 

□  □  Static 


Figure  7:  Predicted  dynamic  and  static  bearing  reaction  forces  versus  shaft  speed,  (a)  M=1  kg,  (b)  M=10  kg 

V.  Conclusions 

The  paper  details  the  unsteady  bulk-flow  transport  equations  for  prediction  of  the  transient  forced  response  of 
turbulent  flow  hydrostatic  bearings,  the  scheme  for  numerical  solution,  and  including  its  integration  to  the  dynamic 
response  of  a  simple  rotor-bearing  system.  The  computational  bulk-flow  model  extends  earlier  analyses  and  enables 
the  current  computational  code  for  ready  use  in  the  prediction  of  transient  rotor-speed  start  ups  in  cryogenic 
turbopumps. 

Predictions  for  the  transient  response  of  a  point  mass  supported  on  one  water  lubricated  six-pocket  hydrostatic 
bearing  are  advanced.  A  schedule  of  a  fast  rotor  speed  start  up  is  specified  along  with  variations  in  feed  supply 
pressure  into  the  bearings  and  increasing  static  load  acting  on  the  rotor.  A  linear  rotordynamic  analysis  predicts  rotor 
unstable  motions  with  ~50%  whirl  frequency  subsynchronous  motions  for  rotor  masses  exceeding  a  critical  value. 
Nonlinear  transient  rotor  response  predictions  are  obtained  for  two  rotor  masses,  one  low  and  the  other  nearly 
equaling  the  critical  magnitude.  For  the  case  with  largest  rotor  mass,  the  predictions  evidence  the  onset  of  unstable 
(subsynchronous)  whirl  motions.  On  the  other  hand,  the  predictions  derived  for  the  low  rotor  mass  show  a  quasi¬ 
static  dynamic  force  response  with  little  effect  from  the  rotor  accelerations  and  no  unstable  motions. 
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(a)  rotor  mass,  M=  1  kg 


FB 


(b)  rotor  mass,  M=  10  kg 


FB 


Figure  8:  External  force  and  predicted  bearing  reaction  force  versus  time,  (a)  M=1  kg,  (b)  M=10  kg 
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Hydrostatic  bearings  for  turbopumps 


Low  cost  primary  power  cryogenic 
turbo-pumps  (TP)  are  compact, 
operate  at  high  speeds,  and  require  of 
externally  pressurized  fluid  film 
bearings  to  support  radial  and  thrust 
loads. 


Thrust  hydrostatic  bearing 
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Radial  hydrostatic 
Bearing 


i  lyunu  u  ii  uo  i  ex  iouioi  ucdiinyb  clldUlc 

smaller  and  lighter  turbopumps  with  no 
DN  life  limitations 

Large  stiffness  (accuracy  of 
positioning)  and  damping  force 
coefficients  allow  for  unshrouded 
impellers  with  increased  TP  efficiency 


Justification 
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Hydrostatic  Bearings  for  Cryogenic  Turbo  Pumps 


Major  Features: 

1  One  piece  titanium  rotor. 

2  Split  hydrostatic  bearings. 

3  Cast  pump  housing  with  integral  crossover  passages. 

4  Radial  inflow  turbine. 

5  Cast  turbine  housing  with  vaneless  inlet  volute. 

6  Filtered  bearing  supply. 

Minck,  A.,  and  Peery,  S.,  1998,  “Design  and  Development  of  an  Advanced 
Liquid  Hydrogen  Turbopump,”  AIAA  paper  No.  98-3681, 34th 
ASAA/ASME/SAE/ASEE  Joint  propulsion  Conference  &  Exhibit 


Advanced  Liquid  Hydrogen  Turbopump 
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Hydrostatic  Bearings  for  pump  applications 


External  pressure  source  forces  fluid  to  flow  between  two 
surfaces,  thus  enabling  their  separation  and  the  ability  to 
support  a  load  without  contact. 

Support  very  large  loads.  The  load  support  is  a 
function  of  the  pressure  drop  across  the  bearing  and 
the  area  of  fluid  pressure  action. 

Load  does  not  depend  on  film  thickness  or  lubricant 
viscosity 

Long  life  (infinite  in  theory)  without  wear  of  surfaces 

Provide  stiffness  and  damping  coefficients  of  very 
large  magnitude.  Excellent  for  exact  positioning  and 
control. 
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Advantages  of  Hydrostatic  Bearings 
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Hydrostatic  Bearings  for  pump  applications 


Require  ancillary  equipment.  Larger  installation 
and  maintenance  costs. 


Need  of  fluid  filtration  equipment.  Loss  of 
performance  with  fluid  contamination. 

High  power  consumption:  pumping  losses. 

Limited  LOAD  CAPACITY  -  f(Psuppiy) 

Potential  to  induce  hydrodynamic  instability  in 
hybrid  mode  operation  (with  rotor  spinning). 


orifice 

i  x 

journal 


Potential  to  show  pneumatic  hammer  instability 
with  compressible  fluids 
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Disadvantages  of  Hydrostatic  Bearings 
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Conventional  hydrostatic  bearing  design 


A 


*  large  pocket  area  (80-90  %  of  total  area) 

*  deep  pocket  depth  (»>  clearance) 

*  large  orifice  discharge  volume 

Applications: 

low  or  null  surface  speed,  low  frequencies, 
nearly  incompressible  fluids  (water  or  mineral  oil) 
produces  very  large  DIRECT  Stiffness. 

Warning:  This  design  should  NOT  be  used  with 
compressible  liquids  or  gases 
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Hydrostatic  Bearings  for  Cryogenic  Turbo  Pumps 


*  small  pocket  area  (10%-25  %  of  total  area) 

*  shallow  pocket  depth 

*  small  or  null  orifice  discharge  volume 

Applications: 

high  surface  speeds,  low  and  high  frequencies, 
compressible  liquids  (LO2,  LH2,  LN2) 


Nearly  inherent 
restrictor  type, 
i.e.  orifice 
coefficient 
regulated  by 
clearance 


+  Angled  injection  against  rotation 
to  reduce  cross-coupled  stiffnesses 
(avoid  hydrodynamic  instability) 
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Cryogenic  fluid  hydrostatic  bearing  design 
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Past  Work  at  TAMU 


HYDROJET®- bearing  model  (1990-1996) 

hydrostatic/hydrodynamic  radial  bearings, 
angled  injection,  roughened  surfaces 

Honeycomb  seals  and  annular  damper  seals 

tilting  and  flexure  pivot  journal  bearings, 
simple  foil  bearings, 

HYDROTHRUST®  -  bearing  model  (1998-2000) 

hydrostatic/hydrodynamic  thrust  bearings 

inner  pressurized  face  seals  with  angular  misalignment 


Bulk-Flow  Codes  include  full  fluid  inertia,  turbulence  flow  and 
thermohydrodynamic  models  for  high-speed,  high-pressure, 
hot/cold  cryogenic  and  process  fluid  operating  conditions. 
Cryogenic  fluids:  LO2,  LH2,  LN2  &  Methane 
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Past  Work  at  TAMU  (Childs  et  al.) 


Hydrojet®  -  radial  hydrostatic  bearings 

Tests  with  water  (1000  psi  max,  25  krpm  max). 

+  20  bearings  x  3  clearances  &  2  pocket  depths,  different  pocket 
shapes,  macro-roughness  (surface  textured)  bearings,  angled 
injection. 

Water  Lomakin  Bearings 
Gas  Honeycomb  seals 

Mineral  Oil  tilting  and  flexure  pivot  journal  bearings 


Hydrothrust  ®-  axial  thrust  hydrostatic  bearings 

NONE  available  in  literature  for  high  speed,  high  pressure  (turbulent 
flows).  Tests  planned  for  Su  2007 


Validation  (benchmarking)  of  bulk-flow  codes 
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Past  Work  at  TAMU  (Childs  et  al.) 


Physical  variable  - 
prediction 

Experimental 
validation  (  +/- 
Accuracy  %) 

Notes 

Bearing  flow  rate  6t  pocket 
pressures 

YES  (2  %) 

Empirical  orifice  coefficient  extracted 

Drag  torque  6t  power  and 
temperature  raise 

YES  (2  %) 

Must  not  reduce  TP  efficiency.  Cryogen 
cannot  exit  too  hot  (vaporization) 

load  capacity  (fluid  film 
forces  and  restoring 
moments), 

YES  (5  %) 

Radial  Load  -  Linear  with  journal 
displacement 

16  complex  impedance 
force  coefficients  due  to 
dynamic  journal  center 
displacements  and  journal 
axis  rotations 

YES 

(5  %,  20%) 

(10%) 

(30  %) 

Lateral  (radial)  motions  only 
stiffness:  Kxx,  KYY  ,  KXY,  KYX, 
damping:  Cxx,  CXY,  CYX, 

inertia:  Mxx,  MYY  MXY,  MYX, 

Whirl  frequency  ratio  for 
lateral  shaft  motions 

YES  (10%) 

WFR  ~  0.50  for  smooth  surface  bearing 
and  radial  fluid  injection 

Warning  for  pneumatic 
hammer  instability 

YES 

Tests  with  water-GN2 

Test  Validation  of  Radial  Hydrostatic  Bearing  Predictions 
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Objective 


Prediction  of  Start-up  response  of  an 
all-fluid  film  bearing  cryogenic  TP  by 


creating  capability  for  modeling  non-linear 

forced  response  of  fluid  film  radial 

bearing,  i.e.  bearing  reaction  forces  (impedance 
models)  as  a  function  of  instantaneous  journal 
position,  velocity  and  acceleration 


Procedure:  Solution  of  unsteady  bulk-flow 
equations  in  fluid  film  bearing  and  integration  to 
(nonlinear)  rotordynamics  model  for  prediction  of 
rotor-bearing  transient  response.  11 
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Unsteady  bulk-flow  equations  in  film  lands 


-Turbulent  flow  with  fluid  inertia  effects 
-Mean  flow  velocities  -  average  across  film  (h) 
-Not  valid  in  flow  zones  with  strong  recirculation 


Continuity: 


Axial 

momentum 
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Unsteady  Bulk  Flow  Equations  in  Pockets 


Pr=Pr+MK, 


QR 


2 (h  +  hR)2  {  2 


■V. 
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Pocket  pressures:  angled  &  radial  injection 
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Unsteady  Bulk-Flow  Analysis  of  Hybrid  Bearings 


Find  bearing  reaction  forces  (Fx,Fy)  as  a  function 
of  rotor  position  and  velocity. 

Other  quantities  of  interest  are:  flow  rate,  drag  torque,  power  loss,  fluid 
temperature  rise.  Implements  perturbation  of  unsteady  flow  field  to  calculate 
(stiffness,  damping  and  inertia)  force  coefficients 

^ap  ~~  ^  { Pp  H a  dz  d6  —  Kap  —  (o  M ap  +ico  Cap  ,  a^=Xj 

A B 

Numerical  method  of  solution: 

Fluid  Flow:  SIMPLEC  control-volume  method,  implicit  time 
response  scheme  with  local  linearization 

Rotordynamics:  Integration  into  in-house  rotordynamics 
program  (FE-component  mode  synthesis) 
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Results  &  Numerical  Analysis  Methods 
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Example:  start  up  water  HB  test  rig 

Fluid:  warm  water 

Max  rotor  speed  =  20  krpm, 

*  | 

Applied  loads  = 
Fn.  of  rotor  speed 


Pocket  depth Ic  ~  7 
Pocket  area  ~  20  % 
to  avoid  water  hammer 


Emulates  typical 
start-up  conditions 
in  a  cryogenic 
turbopump 


/>  Bearing  Loads  FROM 
PUMP  Start  Up 


Pressure  supply  -  RPM2 
(bleed  off  from  pump 
discharge) 


Predict  transient  response 
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Geometry  and  operating  conditions  for  water  HB  rig 
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Rotor-Bearing  System  EOMS 


Numerical  integration  of 

M  ax  =  Wx+Fx,  M  ay=  Wy+Fy 

where 

ax=d(vx)/dt,  ay  =  d(vy)/dt :  accelerations 
vx  =  dX/dt,  vy  =  dY/dt:  velocities 

M:  point  mass  (rigid)  rotor 

Wx,  Wy:  external  loads  including  weight 

Fx,Fy:  bearing  reaction  loads 


Investigate  effect  of 
rotor  mass  on 
transient  response 
and  stability  of  rotor¬ 
bearing  system 
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Example:  Start  up  of  Water  Test  Rig 


Feed  pressure 
into  bearings 
and  external 
loads  on  rotor 
tied  to  rotor 
speed 


Supply  pressure  and  static  load  vs.  rotor  speed 
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Static  Load  Wx  (N)  ^  Supply  Pressure  (bar) 


Example:  Start  up  of  Water  Test  Rig 
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Rotor  speed,  supply  pressure  and  static  load  vs.  time  for  water  HB 
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Rotor  eccentricity  (e/C) 


Example:  Start  up  of  Water  Test  Rig 


Static  &  Dynamic  Response 


-  Transient 

□  □  Static 


Rotor  speed  (krpm) 


Locus  of  journal  center 


-  Transient 

Static 


For  low  rotor  mass:  smooth  and  stable  transient  response,  rotor  reaches 
steady  state  condition.  Quasi-static  analysis  is  sufficiently  accurate  for 

prediction  of  response 


Transient  rotor  position  (eccentricity)  vs.  shaft  speed,  M—  1  If  9 
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Rotor  eccentricity  (e/C) 


Example:  Start  up  of  Water  Test  Rig 


Static  &  Dynamic  Response 


-  Transient 

□  □  Static 


Rotor  speed  (krpm) 


Locus  of  journal  center 


-  Transient 

Static 


For  high  rotor  mass:  smooth  and  stable  transient  response,  but  rotor  begins 
to  whirl  at  max  speed.  Quasi-static  analysis  is  NOT  accurate  for  prediction  of 

system  response 


Transient  rotor  position  (eccentricity)  vs.  shaft  speed,  M—  0  kg 
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Dim  rotor  displacements 


Example:  Start  up  of  Water  Test  Rig 


* 

* 


....  y/c 

e/c  For  low  rotor  mass:  smooth  and  stable  transient 
response,  rotor  reaches  steady  state  condition. 


Transient  rotor  position  (eccentricity)  vs.  time,  VI— 1  kg 
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Dim  rotor  displacements 


Example:  Start  up  of  Water  Test  Rig 


—  e/c  For  high  rotor  mass:  Unsteady  oscillations 

lead  to  rotordynamic  instability  at  max  speed 


Transient  rotor  position  (eccentricity)  vs.  time,  A—  10  kg 
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Conclusions 


Predictions  for  the  start-up  response  of  a  rigid  rotor  supported  on  one 
water  lubricated  hydrostatic  bearing.  Schedule  of  rotor  speed  start  up 
specified  along  with  changes  in  supply  pressure  into  the  bearing  and 
increasing  static  load  acting  on  rotor. 

Linear  rotordynamic  analysis  predicts  rotor  unstable  motions  with  ~50% 
whirl  frequency  subsynchronous  motions  for  rotor  masses  exceeding  a 
critical  value. 

Nonlinear  transient  rotor  response  predictions  obtained  for  two  rotor 
masses,  one  low  and  the  other  nearly  equaling  the  critical  magnitude. 

For  largest  rotor  mass,  predictions  verify  onset  &  persistence  of 
unstable  (subsynchronous)  whirl  motions. 

Predictions  for  low  rotor  mass  show  a  quasi-static  forced  response 
with  little  effect  from  rotor  accelerations  and  showing  no  unstable 
motions 
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Hydrostatic  Bearings  -  Recommendations 


HBs  have  WFR  >  ~  0.50  limiting  their  application  to  ~2x 
critical  speed.  Limiting  speed  condition  can  be  LOWER  if 
fluid  is  compressible  and  pockets  are  too  deep  &  of  large 
area.  fo  reduce  risk  of  hydrodynamic  instability  &  to  increase 

bearing  stability  margin: 


-Texture  bearing  surface 

Proven  with  macro  “rough”  surfaces  such  as 

Knurled,  round  hole  and  tire  truck  pattern  tested  successfully 

-Angled  injection  against  rotation 

Retards  circumferential  flow  swirl,  effectiveness  reduces  at  high  rotor  speeds, 
Can  induce  backward  whirl.  Tested  successfully 

-Introduce  bearing  asymmetry 

Geometrically  build  stiffness  orthotropy  (Kxx  >  Kyy) 

Axial  feed  grooves,  mechanical  preload,  etc.  Tested  &  patented!  24 


Approved  for  Public  Release,  distribution  unlimited. 


Hydrostatic  Bearings  -  Learn  more 


Learn  more: 


http://phn.tamu.edu/TRIBGrou 


Thanks  to  AFRL  &  Northrop  Grumman  for  interest  &  financial  support 
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